Peng-Robinson method